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Abstract 



is studied. By means of the complete Schwarzshild tetrads, we again deal with 
the phase in terms of the type-I and type-II phases, which correspond to the 
momentum and "spin connection" of Dirac particle coupling to the curved 



^ . spacetime respectively. We find that the type-II phase is in the magnitude of 

the square of the gravitational potential, which can be neglected in the weak 
field. For the high energy particle (mass neutrinos), we obtain that the phase 
calculated along the null is equivalent to the half phase along the geodesic in 
the high energy limit. Further, we apply the covariant phase to the thermal 
neutron interference, and obtain the consistent interference phase with that 
exploited in COW experiment. 
PACS number(s): 95.30. Sf, 26.30.+k, 14.60.Pq 
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I. INTRODUCTION 



The neutrino oscillations have been a hot topic in the high energy experimental and 
theoretical phycics recently |jl|,0], in particular, with highly confident atomspheric neutrino 
experiment of Super-Kamiokande to assure the neutrino mass [0]. As a natural extension 
of the theoretical consideration, the description of neutrino oscillation in the fiat spacetime 
should be replaced by that in the curved spacetime if the gravitational background is taken 
into account. In other words, the physics related to the neutrino oscillation in Minkowski 
spacetime with Lorentz invariant will be extended to Riemanian spacetime with general coor- 
dinate transformation. The pioneer theoretical considerations on the gravitationally induced 
neutrino oscillation were proposed by Wudka and later advocated by Ahluwalia [^] who 
for the first time presented the three fiavor neutrino oscillations in the weak field expansion 
scheme. Further some thoughtful idea on the nongeometrical element in the description of 
the gravitational theory, trigered by the gravitationally induced neutrino oscillation clock, 
was also provided |Q. Moreover the violation of the equivalence principle was also employed 
to account for the significant infiuence on the MSW effect P for the solar neutrinos [^H. 
More recently, from the different angle, the gravitational effects on the neutrino oscillations 
have been paid much attention by a number of authors p|,P-p^, but, unfortunately, debates 
and conflicts occur in the understanding of the gravitationally induced neutrino oscillations 

iilMl • 



In this paper, we will discuss the particle interference phase in a gravitational fleld 
in a unifled version, i.e., provide a unifled description of the phases for both high energy 
particle (two mass neutrinos) and low energy particle (COW thermal neutrons). In order to 
develope the treatments by using the null condition to calculate the neutrino relative phase 
and appling the weak fleld condition to calculate COW thermal neutron phase, we employ 
the accurate particle world line (geodesic) to calculate both phase factors, and at last we 
can obtain the correct neutrino relative phase by dividing a factor of 2 and obtain the COW 
neutron interference phase by the convenient approximated condition. 



2 



The energy condition in the gravitational field to account for the relative phase of mass 
neutrinos often makes the confusions and even confiicts [^,|Tl|,|l^] . Apparently, we point out 
that almost all debates related to the neutrino oscillation in curved spacetime originate from 
the inconvenient use of this condition. It is noted that replacing the null by the geodesic to 
calculate the neutrino oscillation phase will produce a factor of 2 error. However this factor 
of 2 will be automatically deleted when considering the two neutrino arrival time differnce 



1^ . Here, we still follow the plane wave treatment for the extremely relativistic neutrinos 



in the framework of the standard treatment otherwise the wave packet treatment will 
be applied for the general case . 



For the reason of simplicity, we confine our treatment in two generation neutrinos (elec- 
tron and muon) and mainly in Schwarzschild geometry with radial propagation and nonradial 
propagation respectively. The purpose of the paper is threefold. First we point out the gen- 
eral treatment of the high (low) energy particle phase in a gravitational field, include the 
type-I and type-II phases. Second we give the complete description of neutrino phase along 
the geodesic, and point out the relation to that by the null in the high energy limit. Third 
we apply the covariant geodesic phase to the COW thermal interference phase, then we can 
find that the geodesic phase can result in the correct results for the interference phases of 
both the two mass neutrinos and the thermal neutron by the standard treatment. 

So the paper is organized as follows. In Sec. II, we discuss Dirac equation and the 
treatments of particle interference phase in the curved spacetime, which includes type-I and 
type-II phases. In Sec. HI and Sec. IV, we calculate the neutrino oscillation phase in 
Schwarzschild spacetime along the geodesic, and compare it with that along the null, in the 
radial and nonradial directions respectively. The application of the unified geodesic phase 
to the COW thermal neutron interference is done in Sec. V. Furthermore, discussions and 
conclusions involved in the high energy neutrino oscillations and low energy thermal neutron 
interference in curved spacetime are given. We set G = h = c = 1 throughout this article. 
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II. DIRAC EQUATION IN A CURVED SPACETIME 



Previous to entering our main point, we stress that the semi classical approximated Dirac 
particle does not follow the geodesic exactly, but the force aroused by the spin and the 
curvature coupling has a little contribution to the geodesic deviation |T3]. So here we take 



the neutrino as a spinless particle to go along the geodesic The gravitational effects on 
the spin incorporated into Dirac equation through the "spin connection" appearing in the 
Dirac equation in curved spacetime [|l6ip!7|, which is constructed by means of the variation 



of the covariant Lagrangian of the spinor field as, 

[7X(^M + rM) + ^]^ = o. (1) 

In this equation and in the rest of this section, greek indices fi, a refer to general covariant 
(Riemanian) coordinates, while the latin indices a, b, c, d refer to locally Lorentz (Minkowski) 
coordinates. The tetrads ha'^ connect these sets of coordinates by 

= Vab /i% h\ , (2) 

the tetrads are supposed to satisfy the following relation 

h%ha'' = 6^'' ■ h\h^ = 5\. (3) 

The explicit expression for T ^ can be written in terms of the Dirac matrices and tetrads (see 

O) 

T, = \Yl'n^]h'^K,,,. (4) 

We must first simplify the Dirac matrix product in the spin connection term. It can be 
shown that 

7"[7^ 7^=] = 2r7<^S^ - 2r7'^^7^ - 2ze'^'^^^757d, (5) 

where rj"'^ is the metric in fiat space and e"'"^'^ is the (fiat space) totally antisymmetric tensor, 
with e"^^^ = +1. With Eq.(^), the contribution from the spin connection is arranged as 
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= ^v^ - ja^75, (6) 
where is the tetrad vector and is the tetrad axial-vector respectively, defined by 

= T'^Xfj,, 0'^l = -^^uXaT"'^'^ (7) 



(8) 



If Vfj, — > 0, then the expression of Eq.(P) is recovered to the form obtained by Cardal and 



Fuller |TT|]. It is instructive to outline the properties of the spin connection. We can see 
that the two terms in the spin connection represent the different aspects of the gravitational 
coupling with the Dirac field. The second term of r.h.s. of Eq.(0), proportional to ^75, has 



the similar form to the weak interaction field [18,19|. Moreover is an axial-vector, which 



represents the physical modification from axial symmetry to the spherical symmetry [21 



From the gravitational interation point of view, shows the rotational gravitational field 
property, or the angular momentum aspect of gravitational field, in other words, it represents 
the gravitomagnetic like interaction. However the first term of Eq.(^) is similar to the 
canonical momentum, then its contribution seems to be the gravitational field momentum. 
In order to group Eq. (P) with terms arising from the matter effects, we can without physical 
consequence arrange Eq.(||) as follows 



(9) 



(10) 
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'2h. 
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In these equations, h = {—gY^'^ = [det(g^,y)]^/^. The expression in Eq.(|^) treats left- and 
right-handed states differently. Proceeding as in the discussion of matter effects, we will 
borrow the neutrino oscillation standard treatment in a gravitational field in |TT[, where 



the gravitational induced effective mass can be calculated from the mass shell condition, 
which is obtained by iterating the Dirac equation 

{P^ + Pg>. + Ag.Vl) {P^ + P^ + A^c'Pl) = m\ (12) 

where is the 4-momentum, and we have not included background matter effects. However, 
the following approximated conditions are usually valid for neutrinos, i.e., » Pcfi, 
P^ » Aof, and P^ » m. 

The complication in calculating the neutrino phase in curved spacetime is related to the 
nature of the neutrino trajectories. In flat spacetime, the neutrino trajectories are straight 
lines. But in the curved spacetime, the geodesic is curved from the global point of view, 
which is more complicated than the situation of the fiat spacetime. Now we will follow the 
factual mass particle world line to cope with the neutrino oscillation problem, but we can 
prove that the neutrino phase along the null is the half of the value along the geodesic in the 
high energy limit (see APPENDIX A). However the relative phase of two mass neutrinos 
should be equivalent to that along the geodesic through dividing a factor of 2 because the 
arrival time difference of two mass neutrinos should be considered ||T^ . 



In fiat spacetime, the phase factor can be written as a conventional manner |TT],[T^,^, 

$ = y" mds = j{Edt - Pdx) = Jri^^P^dx", (13) 

where phase factor $ is also the classical Langrange for the particle motion, and the route 
is along the particle world line (geodesic) determined by the variational principle or Jacobi- 
Hamilton equation P JT5|JT7[| . ds is the interval in flat spacetime with ds"^ = rj^^dx^dx^ and 
the metric rj^^ = diag{+l, —1, —1, —1). In the curved spacetime, however, the neutrino 
phase is calculated along the null pn]JT3[] 

$ = y" mdX = J g^^P^dx'' = J P/^ ^ (14) 

where A is the affine parameter along the null. For the null (photon trajectory) it therefore 
may be convenient to leave the affine parameter A as the variable of integration in Eq.(|l^). 



The tangent vector to the null = dx^/d\, and x^(A) = [x°(A), X"^(A), x^(A), x^(A)]. Now 
we consider the general mass shell condition and ignore the terms of 0{A^), 0{AM), 0{Pq^ 



and 0{Pgij.M). Following ref. [11 , we obtain 



n 



(15) 



From Eq.(|T5|), we find that the total phase includes two types. The first type is con- 
tributed by the phase of the spinless particle w? /2Po, which is conventionally discussed 



ri| , p!3[| and the second type is contributed by the "spin connection". According to the 



definition by Anandan the type-I phase represents the contribution of the particle 
momentum coupled to the spacetime cuvature, and the type-II phase represents the spin 
connection contribution of the particle coupled to the spacetime geometry. 



Further, it is convenient to define a column vector of fiavor amplitudes [|1T|. For example, 
for mixing between z/g and z/,-, 

(^el^(A)) 

KI^(A)), 

Eq.(p!6l) can be written as a differential equation for the null parameter A, 



(16) 



■ dx 



(mV2Po -Ac.n^- Pg, n ^) x, 



(17) 



Eq.(|T7D can be integrated to yield the neutrino flavor evolution. A similar equation was 



obtained in Ref. [11 . 



III. RADIAL NEUTRINO OSCILLATION IN SCHWARZSCHILD SPACETIME 

In this section we will study the phases along the geodesic and along the null in 
Schwarzschild spacetime. In Schwarzshild geometry, the spherically static spacetime can 
be globally represented by the Schwarzshild coordinate system {x^} = {t, r, 6, (p) with the 
line element as 

ds'^ = goodt^ + giidt^ - r'^dO'^ - r"^ sin^edcj)^ , (18) 
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(19) 



where = 2M is Schwarzschild radius and M is the gravitational mass of star. The 
general form of tetrads in Schwarzshild spacetime are given as follows (see APPENDIX B). 
The spacelike components of coordinates share the SO (3) rotational symmetry and timelike 
coordinate shares the reflection symmetry 

/7oo \ 

'~fiis6c(j) r cdccj) —r s6s(, 
'yus9s(f) r c9s(j) r si 
V 7iic^ -r sO 
where 7oo = ^/g^, 7ii = \/~gii, s = sin and c = cos, and the subscript /i is the column 
index. This tetrad deduces the tetrad axial vector and tetrad vector = and v'^ = 
(0,f^,0,0). It is easy to understand that the axial tetrad is cancelled in the case of the 
spherical symmetry, which represents the deviation from the spherical symmetry and occurs 



(20) 



in the Kerr spacetime |21|. The nonzero tetrad vector(see APPENDIX B) 

, 1 - Vdu 



vi = gnv 
1 



-{go 



[Iriy/lfo 



gii 



^i-gn) +^(i-v/^), 



(21) 



so the type-II phase is 



J J Vn 



r — ro + In- 



'l - 



1 - 



(22) 



From Eg . (p2D , if r — r^, the type-II phase will be infinite. Then in the weak field, we 
can conclude that the type-II phase for neutrino in the curved spacetime is proportional to 
the square of the gravitational potential, i.e., (|^)^, which is lower than the contribution 



of the type-I phase that is related to the gravitational potential, i.e., the gravitational red 
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shift [pT] , pr3[ . But, if we consider the strong gravitational field, e.g., near the horizon of the 
black hole 0, the type-II phase might produce some effects on the particle interference. 
However, we stress that the type-II phase of Eq.(|22|) is independent of the neutrino mass, 
namely it is equally contributed to the different neutrinos. So this leads to the cancelling of 
the relative type-II phase between the different mass neutrinos. Therefore we only consider 
the type-I phase of the mass neutrinos in the Schwarzschild geometry by means of the particle 
geodesic equation 

+ ri.-7--7- = o (23) 



ds"^ ds ds 



where 



= \9^'' [dt.9pu + d,gp^ - d^g^,] (24) 

is the Levi-Civita connection. Along radial direction with d6 = d(j) = 0, there are two non- 
trivial independent equations in the four geodesic equations(see Appendix C). The nontrivial 
equation of timelike 

i+^i = (25) 

9oo 

where the dot represents the derivative to ds, for instance, i = The Eq.(|^) gives 

9oot = const., or, E = Pg = mgoot = const. (26) 

Eq.(p6D represents that the covariant energy is the motion constant along the geodesic, by 
which the calculation of the phase is proceeded. It is stressed that it is the covariant energy 
Po (not P°) the constant of motion, then this fact will be important in the later calculation 
of the neutrino phase. Otherwise the ambiguous definition of the neutrino energy in the 
gravitational field will lead to the confusion in understanding the gravitationally induced 
neutrino oscillation. 

The mass shell condition can be written in the Schwarzshild spacetime, which is not 
independent and can be deducted from the geodesic equations 
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goot^ + gur^ + r^4>'^ = 1, 



(27) 



substituting Eq.(|26|) into Eq.(|27|) with <^ = 0, which means the radial motion, thus we obtain 

(28) 



ds 
dr 



-9oo9n 



■Po\2 



-9o 



If m — > 0, then ^ = 0, the null condition is recovered. If r — > oo, the asympotic fiat 
solution goo 1 and —gu ^ 1 is obtained, by which we can acquire the phase of flat 
spacetime. In Schwarzshild spacetime, —googu = 1, Eq.(^) becomes simply 

1 



ds 
dr 



m)^ 9oo 



(29) 



^(^)^. The accurate phase factor in Schwarzshild spacetime along the radial 



where 5^ 

direction is acquired in virtue of the phase factor 

,ds, , f mdr 



^{geod) 



m[—)dr 
dr 



Po\2 



9o 



,mr . I rg 



,rnro^ l^^r_^ 



+ 



, mrg 
^2fc3/2 
, mrg 

.2 



)log 
)log 



rs + 2kr + 2r\/k{k + 



2kro + 2roJk{k + 



To 



m 



0{6^) = 2<!>{nuU), 



(30) 



where k = (^)^ — 1. ^{null) = — ro) is the neutrino phase along the null given by 



11,13 



It is noted that the Schwarzschild coordinates r and ro are not the applicable physical 



distances, and the physical distance connected coordinate is given by dL = y/—giidr, where 
dL is the physical distance. The quantity < goo < ^ leads to goo ^ k 10^^ for the 
electron neutrino. In other words, the gravitational contribution appears in one part of 10^^ 



as shown in Ref. [Q. This means that the Schwarzschild gravitation has little contribution 
to the neutrino oscillation phase, which is same as that concluded in [|r^,|I^]. The above 
phase calculation is a precisely treatment, and we do not use any approxamations, except the 
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final approximated simplification. However, we admit that the Schwarzschild coordinates 
employed to calculate the phase factor is not the real physical distance although there exists 
some complications between them in the general case. 

In order to illustrate the validity of the geodesic phase treatment, we calculate the 
geodesic phase again from the other method. From the equations ( P7D and (P^, the deriva- 
tive between the time-like and the space-like coordinates. 



dt 
dr 



9o 



[1 



(31) 



=^(1 

9oo 



+ 6' 



If the neutrino mass m — 0, then 6 



0, which leads to the null condition ^ 



(-gii) 

9oo 



that is conventionally used in the standard treatment ||rT],|13l. For the electron neutrino, 
the mass m ~ 1 ev and the energy ~ 1 Mev, 6"^ ~ lO^"*^^ lH;i9|) which is really a little 
quantity. If we set 6 = 0, then Eg . (pT]) represents the null condition, which will produce the 



null phase [Tl|JT3[] . But, if this almost negligible quantity is taken into account, then what do 
the things occur to the phase calculation? Let us follow the very routine phase calculation 
as proceeded in the standard treatment of the neutrino phase in a gravitational field, 



$ = 1 g^.P^'dx" = J {gooP°dt + gnP'dr). 



The mass shell condition in Schwarzshild spacetime 



(32) 



(33) 



and the approximated energy momentum relation from Eq. 



-9n)P' = Jgooi-gn)P" 



-gn) 
goo 2P°' 



substituting (0) and (|3TD into (0), we obtain 



^(geod) 



r.dt 

gooP^'-^+giiP' 
dr 



dr 



0^2 



-9ii)9ooP''5 



'9u)9o 



■m 



2Pn 



dr 



(34) 
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-9ii)9o 



m 

2 



\/{-9n)9o 



■ m 



2Po 



dr 



J \/{-9ii)9oo{^)dr 



(35) 



m 



^{r - To) + 0(5^) = 2<!>{null). 



It is not strange that the geodesic phase is double of the result by the null |]TT|. This also 
indicates that the application of the null condition arises the correct result, and the geodesic 



phase of two neutrinos despises the two neutrino arrival time difference [12 



Now we discuss the physical coordinate problem because the Schwarzschild coordinate 
employed in Eq.(|35D is not the proper physical distance, however, which is given by 

(36) 



rr 

L-Lo= / V-9n dr 

J Vn 



rWl roJl 

r V Tn 



hi(^y/r - Ts + ^/r^ - ln(^Vro - r^, + ^/r^ 



+ Ts 



We consider the case of a weak field, then the physical distance is approximately obtained 



T T , ^ 

L - Lo ~ r - To + — In — 

2 r„ 



(37) 



where L(Lo) is the proper physical distance corresponding to the coordinate r {r^). In the 
case of weak field as well as the conditions AL = L — Lo L and Ar = r — ^ r, we 
obtain 



Ar={l-^)AL=[1 + V{L)]AL 



(38) 



where V{L) = is the gravitational potential expressed by the physical distance. So the 
relative phase of the mass neutrinos is 

Am^ 



A<l> 



2Pn 



[l + ViL)]AL 



(39) 



where Am^ = — is the mass square difference. So the oscillation length can be given 
by the usual manner 
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Lose = 27r— ^ [1 - V{L)] . (40) 

Likewise, as stated by Cardall and Fuller [|ll|], does not represent the neutrino energy 
measured by a locally inertial observer at rest at finite radius, but rather the energy of the 
neutrino measured by such an observer at rest at infinity. It is generally not possible to ex- 
tract a separate "gravitational phase" from the expression (|39|) ; nevertheless, it is clear that 
gravity has an effect on the oscillations of the radially propagating neutrinos, corresponding 
to the gravitational redshift. In the weak field limit one could define a "gravitational phase," 
however, the energy definition and coordinate condition should be apparently declared. 



IV. NONRADIAL PROPAGATION IN SCHWARZSHILD SPACETIME 



In this section we contrast the radially propagating neutrinos with the nonradically 
propagating neutrinos in order to demonstrate how the gravity affects vacuum neutrino 
oscillations. As a further application, the general covariant phase can be applied to solve 
the nonrelativistic thermal neutron interference. For the general nonradial propagation of 
the mass particle, the knowledge of the geodesic should be exploited p2[ (see APPENDIX 
C). In the 9 = n/2 plane. The accurate phase along the azimuthal geodesic can be calculated 
by 



^"'{geod) = J mds = J m{—)dr 



(41) 



mdr 



Po\2 



+ 



,mr 



lk + - + - 



^k + - + - 



)log 
)log 



/ r P 

Ts + 2kr + 2r]lk{k + + ^] 



2kro + 2rJk{k + -+^-^] 



where / is the integral of motion r^^ = J/m = I = const, and J is the angular momentum of 
the particle. If / = 0, the radial geodesic phase of Eq.(|35|) will be recovered. Equivalently, 
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the nonradial phase can be calculated in terms of the azimuthal coordinate 



ds , 



^"■{geod) = / mds = / m{—)d(j) = / m(0) ^d(j) 



mr 



I 



(42) 



then the particle trajectory can be determined by the trajactory equation p2 



(43) 



where e is the ellipse eccentricity, and we set the initial perihelion location at zero and 
neglect the perihelion precession term. Substituting Eq.(^) into Eq.(^), we have, becomes 



^a^geod) = (^)/(0), 



dx 



+ 



(1 + ecosa;)^ 
2ArcTanh[ 
e sin 6 



(44) 
(45) 



e- l)Tan(0/2) j^^2 _ ^^,-3/2 



(e2 - 1)(1 + ecos0) 
For the low energy particle, the eccentricity of the trajectory can be given 



1 + 



^ m{rsm/2y 



(46) 



2E/m = (f^ -), J = mvr. 



In the case of COW experiment, the scale of the instrument is about 5 centimeters (shown 
in FIG.l.), which is much less than the radius of the earth, and so the transverse azimuthal 
angle is very small = 50. Therefore, Taylor expansion gives 

5(t> 6(j) 



m) 



:i + e) 



:i + |l-2^;2i-|)2- 



(47) 



V. THE COVARIANT TREATMENT OF THE THERMAL NEUTRON 

INTERFERENCE 



Now we use the covariant phase to deal with the thermal neutron interference on the 
earth (COW experiment), which is satisfactorily described by the COW interference phase 
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obtained through solving Schoedinger equation induced by the Newtonian gravitational po- 



tential or through the weak field approximation method Our purpose of using the 
covariant method to cope with the low energy thermal neutron interference is to find how 
the spacetime global structure infiuences on the COW interference. Unlike the covariant 
treatment, the Newtonian gravitational potential induced COW interference phase does not 
include the information of the global structure of the curved spacetime, and nonetheless the 
weak field approximation method induced COW interference phase may also neglect 
some characteristics of the curved spacetime. The second motivation of using the covariant 
method to deal with the thermal neutron interference aims at indicating that the covariant 
phase is not only useful for the high energy mass neutrinos but also successfully apphcable 
to the low energy thermal neutron interference. The particle wave behaviours incorporated 
into the curved spacetime will refiect some quantum properties of the particle in Einstein's 
general relativity, which has not yet been well explored until now. Moreover the successfully 
application of the covariant phase to the COW experiement can help us to wipe out the sus- 
pect of the universality of the covariant treatment of the phase on the particle interference 
in any energy bands. 

In FIG.l., the scale of the instrument of COW experiment is some centimeters, which is 
much less than the scale of the gravitational source, the radius of the earth. For the arm 
plane of the COW instrument perpendicular to the surface of the earth, the thermal neutron 
motions are divided into two beams, the one is radial motion along AB (CD) and the other is 
nonradial along AC (BD). We calculate the phase difference along the different routes ABD 
and ACD. Taking the earth as an exact sphere, the Schwarzschild metric induced radial 
velocity is 



- = ^ = ^{E/mY - g,, = Jk^ , (48) 
m as ^ V r 



so the momentum difference from A to B 



[Pi - PDim' = ^ - (49) 
= 2gH + 0{5s5l), 
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where = -^^ ~ 10 and = ~ 10 ^, g = jkr is the gravitational acceleration on 



the surface of the earth and Rca is the earth radius. 



For the nonradial motion of the nonrelativistic thermal neutron beams through AC and 
BD, with the velocity condition 2f^^ > 1 (in the COW experiement, ~ 10^^° |25[ and 
10^^^), from Eq.(^), we have 



— ^ n~j 

r 



m) - - O'^^, (50) 

where vg = and the velocity approximation condition now becomes ^ > 1, which 
yields the transverse phase at radius r, from Eg. (PI) and Eg. (pO|) 

*"(r) = (^)m) = (51) 

where P = mv is the transverse momentum, and the Schwarzschild coordinate r should 



be replaced by the physical distance if applicable, i.e., dL = y/—gudr, where L notes the 
physical distance. In the case of the weak field, the relation between the coordinate difference 
and the physical distance is Ar = (1 + ^)AL. Thus Eg. (^T]) becomes 

*W = ^(1 + ^), (52) 

where H is the arm length of the the COW interference instrument (FIG.l), and we set the 
sguare route, i.e., H=AB=BD=AC=CD in FIG.l. Then the interference phase difference of 
two thermal neutron beams becomes 

A^COW) = <!>{ABD) - <!>{ACD) (53) 
= <l>(i?e + ^)-<f(^e) 

(P) ^ ^271 

where A„ = 27r/ (P) is the thermal neutron wave length with the momentum (P) = (P^ + 
Pb)/2 ~ Pa ~ Pb- This is just the anticipated result in COW experiment [p5|] 

However, we stress two facts here, from the point of view of the covariant treatment of 
Dirac particle. Firstly, the COW thermal neutron phase obtained is a type-I phase because 
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the type-II phase is the square of the gravitational redshift, which can be neglected in the 
weak field condition. Secondly, if the thermal velocity satisfies v < vq, unlike the COW 
phase, the interference will follow a different regulation related to the thermal neutron wave 
length. The above condition cannot be obtained from the weak field expansion method to 



deal with the COW neutron interference PU|. In other words, we can construct the critical 



wave length for COW thermal neutron as Xn = therefore the critical condition for the 

° ^ mvQ ' 

COW interference phase needs the wave length of the thermal neutron must be shorter than 
this critical wave length , i.e.. An < Xg- 

VI. CONCLUSIONS AND DISCUSSIONS 

On the basis of the particle interference phase along the geodesic, the following conclu- 
sions are obtained and demonstrated. 

(1) On the phase of the particle propagating in Schwarzschild spacetime, we exploit the 
geodesic route to calculate it, and our results for the neutrino oscillation relative phase is in 



agreement with that in Refs. p^Jl3[1 if the geodesic phase is divided by a factor of 2, which 
originates from the consideration of the simultaneous arrival of two mass neutrinos for the 
interference However our results are different from that in ref. [§, where the weak 

field expansion is applied. We find that two factors infiuence on our understanding of the 
gravitationally induced neutrino oscillation, i.e., energy definition and physical distance. The 
geodesic phase will produce a factor of 2 error contrasting to the null phase, which is also paid 



attention in the fiat spacetime p3| , and then the energy definition needs the integral constant 
of motion in the gravitational field. Moreover, the physical distance should be transformed 
from the Schwarzschild coordinate system. Otherwise any ambiguous definitions of above 
two conditions will result in confusion or confiict because the mass neutrino phase is very 
sensitive to these conditions. 

(2) Two types of phases are taken into account in the article. Type-I phase is contributed 
by the curved spacetime, which represents the coupling of the momentum of the particle to 



17 



the curved spacetime. Type-II phase is contributed by the couphng of the "spin connection" 
of the particle to the curved spacetime. In the case of Schwarzschild weak field, type-II phase 
is proportional to the square of the effect of gravitational red shift. Namely, type-I phase 
is the first order effect of the gravitational potential and type-II phase is the second order 
effect of the gravitational potential. So we can neglect the type-II phase when dealing with 
the physical problems in the solar system, such as atomspheric neutrino, solar neutrino as 
well as the thermal neutron interference in COW experiment. 

(3) We find that, if the thermal neutron velocity exceeds over a critical velocity, then 
COW interference phase will appear. But, if the thermal neutron velocity is not higher 
than this velocity, the COW phase will be broken. This critical velocity may be the critical 



condition for the validation of the COW interference 25 , however this conclusion cannot be 



derived from the weak field method to deduce the COW interference phase [^. From the 
quantum mechanics language, there exists a critical wave length, and the COW interference 
needs that the thermal neutron wave length must be shorter than this critical wave length, 
An < with Xg = If this conclusion is valid, we could suggest the low velocity ther- 
mal neutron interference experiement to inspect the proposed effect. Moreover the further 
experiement not only detects the Newtonian gravitational potential induced interference ef- 
fect, but also becomes a sensor of the curved spacetime, to inspect the general relativity 
induced COW interference condition. 

In summary, we can speculate that the type-II phase should play an important role in the 
strong rotational gravitational field, especially near Kerr black hole, where the high order 
terms of gravitational potential will effect. 



APPENDIX A 

The velocity of an extremely relativistic neutrino is nearly the speed of light. Despite of 
this, the propagation difference between a massive neutrino and a photon can have important 
consequences. In the standard treatment of the neutrino oscillation, the neutrino is usually 
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supposed to travel along null-lines [p],p|, [l8| , p!9| , and almost no attention has been paid to this 



small difference. Although seemingly irrelevant, this tiny deviation becomes important for 
the understanding of the neutrino oscillation. Motivated by this argument, we will compare 
the neutrino phase when calculated along the geodesic and along the null-line. With this, 
we will be able to verify the factor of 2 error when the null is replaced by the geodesic. This 
study can be shown to remain valid in the case of flat spacetime. 

Let and n'^ be the tangent vectors to the geodesic and to the null-line, respectively, 
their difference e'^ being a small quantity for the case of an extremely relativistic neutrino. 
Here, we suppose that the two neutrinos, the massless and massive, start their journey at 
the same initial spacetime position A, and their space routes are almost the same. But, 
there will be an arrival time-difference at the detector position B. This means that their 4- 
dimensional spacetime trajectories are not the same, and consequently the tangent vectors 
will present a small difference. Thus, we have 

n'' = + e", or P'' = p'' + me" , (54) 

O o 

where = mn^ {p^ = mn^) is the 4-momentum along the geodesic (nuU-hne) with 

"=^ = ^ 

and 



(56) 



ds 

In these expressions, A and s are respectively affine parameters along the null and the 
geodesic lines. These two tangent vectors satisfy the mass shell relations of the geodesic and 
the null-line: 

g^^n^n" = 1 (57) 

and 

g^.^n" = . (58) 
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Now, substituting (|^) into (|57D, we obtain 

g,u{n^ + e^^){n'^ + e'') = l, (59) 

or, by using (pS]), 



2g^,n^^e^ + 0{e^) = 1 . (60) 

We can estimate the order of {n^} and {n'^} by noting that ne ~ 1/2, which imphes that 
e ~ ~ ^, where E = P° {i = 1, 2, 3) for a relativistic neutrino. 

The neutrino phase induced by the null condition, as in the standard treatment, comes 
from the 4-momentum defined along the geodesic line, and the tangent vector {n^} to 



the null-line |[TT|. We notice that, if the 4-momentum P'^ defined along the null-line was 
instead used to compute the null phase, we would obtain zero because of the null condition. 
Therefore, the phase along the geodesic line (geodesic phase) and the phase along the null- 



line (null phase) can be written respectively as |ri|, p^ , pO| ,p!^ 

$(geod) = j mds = J g^^P^rfds , (61) 

and 

$(null) = j g^yP^'n^ds . (62) 

Therefore, the difference between the geodesic phase and the null phase, by using Eq.(|60|), 
is 



$(geod) - <l>(null) = j gf.uP^in" - n'')ds 

= J g^.P^t^ds = j g^.P^e'^ds + 0{e 
= ^Jmds + 0{e^) = ^$(geod) + 0( 



that is 



$(geod) = 2<l>(null) + 0{e^) . (63) 
20 



This conclusion, valid for a general curved spacetime, is similar to that found in in a 



Schwarzschild [12| spacetime. Concerning the Schwarzschild spacetime, Bhattacharya et 



al |T2| have the following argument for the factor of 2. As the neutrino energy is fixed, 
but the masses are different, if a interference is to be observed at the same final spacetime 
point B(rB,tB), the relevant components of the wave function could not both have started 
at the same initial spacetime point A(ryi,tyi) in the semiclassical approximation. Instead, 
the lighter mass (hence faster moving) component must either have started at the same time 
from a spatial location r < r^, or (what is equivalent) started from the same location ta at a 
later time Ia + At. Hence, there is already an initial phase difference between the two mass 
components due to this time gap, even before the transport from ta to which leads to the 



phase $(null), i.e., the additional initial phase difference may be taken into account [12 



APPENDIX B 



Tetrad in the spherical, static and isotropic coordinate system (X*^ = t,X^,X^,X^) 



^'■^''o = yC{p) and h'^a = y D{p)6^ p6|, the corresponding line element (to avoid confusion, 
Latin indices are enclosed in parentheses to express Lorentz coordinates) 

ds'^ = g^^^dX^dX" (64) 
= C{p)dt^ - D{p){dp^ + p^dQ^) 

rffi^ = de"^ + sin ed(f)'^, X^ = p sin 9 cos (65) 
X^ = p sin 6 sin 0, X^ = p cos 9 (66) 

comparing the line element in the Schwarzschild form, we get 



C{p)=goo, ^lD{p)p = r, ^ = ^1^^ (67) 



using the general coordinate transformation 



h\ = y^h\ (68) 
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(69) 



where {X^} and {X '^} are the isotropic and Schwarzschild coordinates respectively. 
Through using the above transformation and relationship, then we obtain the tetrad in the 
Schwarzschild coordinate system. 

floo \ 

•yusOccp r c6c(j) —r sOscp 
7iis6's0 r cOsip r s6c(j) 
V -fuc9 -r s9 / 
where joo = ^j9oo-, 7ii = \/—9\\-, s = sin and c = cos, and the subscript fj, is the column 
index. The contravariant tetrad can be constructed by /i^^ = g^^h^uVabi 

(Too \ 

l\\s9c<\) c9c(f) — (r s9)~^s(f) 
Jiis9s(f) r^^ c9s(f) (r s9)^^c(l) 
V Tic9 -r--^ s9 ) 

One can proof that the orthagonal relation between h°-^ and h'^^ are satisfied. Eqs. (|69|) and 
(ffOl) can construct the tetrad vector and the tetrad axial-vector 



(70) 



r°oi = [/n^^', 1^32 = 1^23 = tan^ 



r\i = [/nv/^', 
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\r-9 
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r'33 = r^22(sin^)2, 1^33 = -sin ^ cos ^ 



J- 21 — J- 31 — ~) J- 12 — -L 31 — 

d^J-gll 



V-5'11 



(ir 



(71) 



h = r sin 6*, 



^12 
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-T221 = (1 - V^i)/r, 



(72) 



a'^ = 



(73) 



1 - v^^^ 



1 ,r 



s\2 



2r r 
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vi = gnv 
1 



1 



(74) 



911 



^{-9ii) +^(l-v^) 



APPENDIX C 



In the Schwarzschild geometry, the components of the geodesic equation are |22 



i+^i = (75) 

9oo 

r + l9'ooi-9'')iir+ (76) 
ln\[^^ - rg^^df + sin^ e{^f] = 

e + '^Or - sin 9 cos e{^f = (77) 

+ -0r + 2 cot 66^ = (78) 

r 

where tlie dots indicate derivatives with respect to s, such as t = ^ and r = ^. We will 
assume that the orbit is in the plane 6* = |, then Eqs.(|7|) and ( [78D can now be integrated 
directly, giving two integrals of the motion 

• E , , 

9ooi = — = const, (79) 
m 

E = Po, = — = I = const, J = — P3 
m 

These two constants are proportional to the energy and the angular momentum, respectively. 
Precisely, these quantities are constants of the motion and represent the energy and the 
angular momentum per unit mass. 
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FIGURES 




FIG. 1. Diagram illustration of the thermal neutron interference route. We assume the arms 
AB, BD, AC and CD to be same, i.e., AB = BD = AC = CD = H. OA and OC represent the 
radius of the earth and the angle between OA and OC is d(f). 



26 



